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Abstract 

We investigate the problem of secure communication over parallel relay channel in the presence 
of a passive eavesdropper. We consider a four terminal relay-eavesdropper channel which consists of 
multiple relay-eavesdropper channels as subchannels. For the discrete memoryless model, we establish 
outer and inner bounds on the rate-equivocation region. The inner bound allows mode selection at 
the relay. For each subchannel, secure transmission is obtained through one of two coding schemes at 
the relay: decoding- and-forwarding the source message or confusing the eavesdropper through noise 
injection. For the Gaussian memoryless channel, we establish lower and upper bounds on the perfect 
secrecy rate. Furthermore, we study a special case in which the relay does not hear the source and 
show that under certain conditions the lower and upper bounds coincide. The results established for the 
parallel Gaussian relay-eavesdropper channel are then applied to study the fading relay-eavesdropper 
channel. Analytical results are illustrated through some numerical examples. 
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I. Introduction 

In conventional point-to-point wired networks, security is facilitated by secret key sharing 
between relevant parties based on some common cryptographic algorithm. The premise is that 
only legitimate users have access to the encrypted messages and extraneous users (adversaries) 
are unable to access any useful information. The wireless channel is characterized by its inherit 
randomness and broadcast nature. Physical layer security exploits the basic attributes of the 
wireless channel for instance, difference of the fading gains between the legitimate channel 
(source to the legitimate receiver) and the channel to the adversary, to transmit information 
securely to the legitimate receiver. Thus, it eradicates the need of secret key sharing. 

The wiretap channel introduced by Wyner is a basic information-theoretic model which 
incorporates physical layer attributes of the channel to transmit information securely [1]. Wyner's 
basic model consists of a source, a legitimate receiver and an eavesdropper (wiretapper) under 
noisy channel conditions. Secrecy capacity is established when the eavesdropper channel (the 
channel from the source to the eavesdropper) is a degraded version of the main channel (the 
channel from the source to the legitimate receiver). The discrete memoryless (DM) channel 
studied by Wyner is further extended to study some other channels for which secrecy capacity 
is established, i.e., broadcast channels (BC) 0, J3), multi-antenna channels flU-O, multiple 
access channels BH-flU, fading channels ifTOl . [fTT| etc. The idea of cooperation between users 
in context of security was introduced by Ifl2ll . The intuition is that, when the main channel 
is more noisy than the channel to the eavesdropper, cooperation between users is utilized to 
achieve positive secrecy capacity. Secrecy is achieved by using the relay as a trusted node that 
facilitates the information decoding at the destination while concealing the information from the 
eavesdropper. A special case in which there is a physically degraded relay-eavesdropper channel 
was studied in [fT3l . The case in which the relay does not acts as a trusted node is studied in 

ma, ma. 

In this paper, we study a parallel relay-eavesdropper channel. A parallel relay-eavesdropper 
channel is a generalization of the setup in |[T2l|. in which each of the source-to-relay (S- 
R), source-to-destination (S-D), source-to-eavesdropper (S-E), relay-to-destination (R-D) and 
relay-to-eavesdropper (R-E) link is composed of several parallel channels as subchannels. The 
eavesdropper is passive in the sense that it just listens to the transmitted information without 
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modifying it. We only focus on the perfect secrecy rate, i.e., the maximum achievable rate at 
which information is reliably sent to the legitimate receiver, and the eavesdropper is unable to 
decode it. 

The parallel relay-eavesdropper channel considered in this paper relates to some of the channels 
studied previously. Compared to the parallel relay channel studied in |[T6ll . the parallel relay- 
eavesdropper channel requires an additional secrecy constraint. The parallel relay-eavesdropper 
channel without relay simplifies to a number of channels discussed previously. For example, 
the parallel wiretap channel studied in ifTTl . the parallel broadcast channel with confidential 
messages (BCC) and no common message studied in [[3J . 

Contributions. The main contributions of this paper are summarized as follows. For the 
discrete memoryless case, we establish inner and outer bounds on the rate-equivocation region 
for the parallel relay-eavesdropper channel. The inner bound is obtained through a coding 
scheme in which, for each subchannel, the relay operates either in decode-and-forward (DF) 
or in noise forwarding (NF) mode. We note that establishing our outer bound for DM case is 
not straightforward and it does not follow directly from the single-letter outer bound for the 
relay-eavesdropper channel developed in |[T2l Theorem 1]. Therefore a converse is needed. The 
converse includes a re-definition of the involved auxiliary random variables, a technique much 
similar to the one used before in the context of secure transmission over broadcast channels 0. 

For the Gaussian memoryless model, we establish lower and upper bounds on the perfect 
secrecy rate. The lower bound established for the Gaussian model follows directly from the 
DM case. However, we note that establishing a computable upper bound on the secrecy rate 
for the Gaussian model is non-trivial, and it does not follow directly from the DM case. In 
part, this is because the upper bound established for the DM case involves auxiliary random 
variables, the optimal choice of which is difficult to obtain. In this work, we develop a new upper 
bound on the secrecy rate for the parallel Gaussian relay-eavesdropper channel. Our converse 
proof uses elements from converse techniques developed in Q, [0 in context of multi-antenna 
wiretap channel; and in a sense, can be viewed as an extension of these results to the parallel 
relay-eavesdropper channel. This upper bound is especially useful when the multiple access part 
of the channel is the bottleneck. We show that, in contrast to upper bounding techniques for 
our model that can be obtained straightforwardly by applying recent results on multi-antenna 



TO APPEAR IN IEEE TRANSACTIONS ON INFORMATION FORENSICS AND SECURITY 



3 



wiretap channels [@]|-[[6]|, our upper bound shows some degree of separability for the different 
subchannels. 

We also study a special case in which the relay does not hear the source, for example due to 
very noisy source-to-relay links. In this case we show that under some specific conditions noise- 
forwarding on all links achieves the secrecy capacity. The converse proof follows from a new 
genie-aided upper bound that assumes full cooperation between the relay and the destination, and 
a constrained eavesdropper. The eavesdropper is constrained in the sense that it has to treat the 
relay's transmission as unknown noise for all subchannels, an idea used previously in context of a 
class of classic relay-eavesdropper channel with orthogonal components [18]. These assumptions 
turn the parallel Gaussian relay-eavesdropper channel into a parallel Gaussian wiretap channel, 
the secrecy capacity of which is established in 0, ifTTl . 

Furthermore, we study an application of the results established for the parallel Gaussian relay- 
eavesdropper channel to the fading relay-eavesdropper channel. We assume that perfect non- 
causal channel state information (CSI) is available at all nodes. The fading relay-eavesdropper 
channel is a special case of the parallel Gaussian relay-eavesdropper channel in which each 
realization of a fading state corresponds to one subchannel. We illustrate our results through 
some numerical examples. 

The rest of the paper is organized as follows. In section II, we establish outer and inner 
bounds on the rate-equivocation region for the DM channel. In section III, we establish lower 
and upper bounds on the perfect secrecy rate for the Gaussian model, and consider a special 
case in which under some specific conditions secrecy capacity is achieved. In section IV, we 
present an application of the results established in section III to the fading model. We illustrate 
these results with some numerical examples in section V. Section VI concludes the paper by 
summarizing its contribution. 

Notations. In this paper, the notation Xnn is used as a shorthand for (X 1 ,X 2 , . . . ,X L ), the 
notation XJ\ L j is used as a shorthand for (X™, X^ - , . . . , X£) where for / = 1, . . . , L, X™ := 
(Xji,Xj2, • • • , Xin), the notation Xu lu is used as a shorthand for (X^j, X 2j j, . . . , X^j), the 
notation Xm n is used as a shorthand for X n x X 12 . . . x X 1L , E{.} denotes the expectation 
operator, \X\ denotes the cardinality of set X, L denotes the number of subchannels, the boldface 
letter X denotes the covariance matrix. We denote the entropy of a discrete and continuous 
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Fig. 1. The parallel relay-eavesdropper channel. 



random variable X by H(X) and ft.(A') respectively. We define the functions C(x) = \ log 2 (l+x) 
and [x} + = max{0, x}. Throughout the paper the logarithm function is taken to the base 2. 

II. Discrete memoryless channel 

In this section, we establish outer and inner bounds on the rate-equivocation region for the 
discrete memoryless parallel relay-eavesdropper channel. 

A. Channel Model 

Definition 1: The parallel relay-eavesdropper channel consists of four nodes, a source, a 
relay, a destination (legitimate receiver) and a passive eavesdropper. The communication takes 
place over L subchannels. Fig. 1 represents the studied model. The source wishes to send 
confidential messages to the destination, with the help of the relay to conceal them from passive 
eavesdropper. The source encodes the confidential message W to X™ 2 , . . . , X™ L ) codewords 
and broadcasts it over L subchannels to the relay and the destination. The relay helps to reduce 
the uncertainty about the confidential message at the destination by re-encoding whatever it 
has received from the source and transmits {X^.X^ ■ ■ ■ ^21) codewords to the destination. 
The outputs at the relay and destination are given by (Y{[, Y^, . . . , Y^ L ) and (Y", ■ ■ ■ , Y£) 
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respectively. The passive eavesdropper overhears to the source and relay transmission over the 
multiple-access link, which is denoted by (F 2 ", Y^, . . . , YJfo). 

More precisely, the parallel relay-eavesdropper channel consists of <^2[i,l] as finite input 

alphabets and 3^1,-L] a s finite output alphabets. Since the channel is memoryless, 

the transition probability distribution is given by 

L n 

JJ IJp(?/M) !/iM> fa,* I a;iM> x 2m) (1) 
z=i i=i 

where £i M e X lh x 2 i,i G Af M ,j/i,i G 3>j,2/u,i G y u and y a ,i e 3>2i, for / = 1, . . . , L and 
z = 1, ... ,n. The symbols a^j and x 2 i are the source and relay inputs on subchannel I, and 
Uih Vh Vn are the channel outputs at the relay, destination and eavesdropper for the Z-th subchannel 
respectively. 

Definition 2: The source sends a message W E W = {l,...,2 nR } using a (2 nR ,n) code 
consisting of 

• a stochastic encoder at the source that maps W — > L i , 

• a relay encoder that maps /iC^iV n) ~~ ^ ^2[i,L],i f° r 1 < i < ^, 

• a decoding function #(.), that maps the received codewords from the source and relay node 
to get an estimate of the confidential message, g(Yft L S) — > W. 

Definition 3: The average error probability of a (2 nR , n) code is defined as 



Due to the openness of the wireless medium, the eavesdropper listens for free to what the source 
and relay transmit. It then tries to guess the information being transmitted. The equivocation rate 
per channel use is defined as R e = if (W^Y^ L i)/n. Perfect secrecy for the channel is obtained 
when the eavesdropper gets no information about the confidential message W from L y That 
is, the equivocation rate is equal to the unconditional source entropy. 

Definition 4 ( fiTfl): A rate-equivocation pair (R, R e ) is achievable for the parallel relay- 
eavesdropper channel, if for any e > there exists a sequence of codes (2 nR , n) such that 




(2) 
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for any n > n(e) 

n 

n 

P: < e. (3) 

B. Outer Bound 

The following theorem provides an outer bound on the rate-equivocation region for the parallel 
relay-eavesdropper channel. 

Theorem 1: For a parallel relay-eavesdropper channel with L subchannels, and for any achiev- 
able rate-equivocation pair (R,R e ), there exists a set of random variables Ui — > (Vu,V 2 i) — > 
(X ll} X 2l ) (Y h Y lh Y 2l ), I = 1, . . . , L, such that (R, R e ) satisfies 

L L 

R < minl^/^,^;^),^/^;^,^ | V 2l )] 
i=i i=i 

Re — Ft 

L 

R e < min { I(Vu, V 2l ; Y t \ U t ) - I(V lh V 2l] Y 2l | 
i=i 

L 

J2l(Vu;Y h Y u | Va,Ui) - I(V lh V 2l ;Y 2l \ Ui)}. (4) 
i=i 

Proof: The proof of Theorem 1 is given in Appendix HI □ 

Remark 1: The outer bound in Theorem 1 does not follow directly from the single-letter 
outer bound on the rate-equivocation region established for the relay-eavesdropper channel |[T2l 
Theorem 1]. Therefore a converse is required, in which we need to re-define the involved auxiliary 
random variables. The technique used to re-define the auxiliary random variables has some 
connection with the one used before in the context of secure transmission over broadcast channels 
0. 

Remark 2: The region © reduces to the rate-equivocation region developed for the relay- 
eavesdropper channel |[T2l Theorem 1] by setting L := 1 in ©. 

Remark 3: The equivocation rate in Theorem 1 reduces to the secrecy capacity of the parallel 
wiretap channel established in j3l Corollary 1] by removing the relay, i.e., by setting Yu = X 2 i = 
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V 2[ = <p. The resulting term ^4=1 HYu', Yi I Ui) — I(Vu) Y 2i | E/j) is maximized by U[ :=constant, 
for I — 1, . . . , L. 

C. Achievable Rate-Equivocation Region 

In this subsection we establish an achievable rate-equivocation region for the parallel relay- 
eavesdropper channel. The achievable region is established by the combination of two different 
coding schemes, namely decode-and-forward and noise forwarding. In DF scheme, for each 
message source associates a number of confusion codewords, the relay after receiving the 
source codewords, decode it and re-transmits it towards the legitimate receiver and eavesdropper 
(see |[T2l Theorem 2] for details). In the NF scheme the relay does not decode the source 
codewords, but transmits confusion codewords independent from the source codewords, towards 
the legitimate receiver and the eavesdropper (see lfT2l Theorem 3] for details). 

Theorem 2: For a parallel relay-eavesdropper channel with L subchannels, the rate pairs in 
the closure of the convex hull of all (R, R e ) satisfying 

R < minl^/^^s^^l^^^'^l^'^^l + S 1 ^^!^ 

leA leA leA c 

Re ^ R 

R e < min { HV11, V 2 f, nU{) - I(V lh V 2l ; Y 2l \U{), ^ I(V U ; Y xl \V 2h U{) - I(V U , V 2l ; Y 2l \U^ 
leA leA 

+ I(V U \ nV 2l ) + min { I{V2U Yi), £ Y 2l \V u )} - min { ^ I(V 2l ; Y{), 
leA c leA c leA c leA c 

J2l(V 2l ;Y 2l )} -J2l(V u ;Y 2l \V 2l ), (5) 

leA c leA c 

for some distribution p(m , vu ,v 2t ,xu, x 2i , yi , yu , y 2 i) = p(ui)p(vu , v 2 i\ui)p(x u , x 2 i\vu , v 2t )p(yi , yu , y 2i \xu , x 2t ) 
for I £ A and p(vu,v 2 i,xu,x 2 i,yi,yu,y 2 i) = p{vu)p(v2i)p(xu\vu)p(x2i\v2i)p(yi,yii, y 2 i\x u , x 2i ) for 
I £ A c , are achievable. 

Outline of Proof: 

The region in Theorem [2] is obtained through a coding scheme which combines appropriately 
DF and NF schemes. In the statement of Theorem [2l sets A and A c represent the subchannels 
for which relay operates in DF and NF mode, respectively. The rates for the DF scheme can be 
obtained readily by setting U := U^^Vi := V 2 := V 2 [i,\a\], Y := Y [UAl] , Y x := Y l[hlAl] 
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and Y 2 := IWuiis f° r I & A in [fT2l Theorem 2]. Similarly the rates for NF scheme can be 
readily obtained by setting V x := Vi[i,|./4*|], V 2 := V 2[lt \ A ^Y := Y [UA c lh Y 1 := ^i[i,|^|] and 
Y 2 ■= Y 2[i,\A c \l> for 1 e A ° in lEl Theorem 3]. □ 

Remark 4: For a parallel relay-eavesdropper channel in which all subchannels are degraded^, 
i.e., 

v{Vh Vih V2i\x\u x 2 i) = p(yu\x lh x 2 i)p(yi\y u , x 2 i)p{y 2i \y h y u , x lh x 2i ), 
I = 1, . . . , L, the perfect secrecy capacity is given by 

L 

C s =maxmm{^[/(^,y 2/ ;H | U t ) - I(V lh V 2l ;Y 2l \ Ut)] + , 
i=i 

L 

J2[I(Vn> Yu | V 2l , Ui) - I(V lh V 2l ; Y 2l \ + } (6) 
i=i 

where the maximization is over U t ->■ (V lh V 2t ) -)■ (X^, X 2/ ) -> (F/, Y lh Y 2i ), for Z = 1, . . . , L. 

Proof: The achievability follows from Theorem 2 by setting A c := 0. The converse follows 
along the lines of Theorem [TJ and is omitted for brevity. □ 

III. Gaussian Channel 

In this section we study a parallel Gaussian relay-eavesdropper channel. Fig. [2] depicts the 
studied model. We only focus on the perfectly secure achievable rates, i.e., (R,R e ) = (R,R). 

A. Channel Model 

For a parallel Gaussian relay-eavesdropper channel, the received signals at the relay, destination 
and eavesdropper are given by 





— X\ij - 






Y lti 


= x u>i - 






Y 2lji 


= X u ,i - 


- \/p~2lX 2 i ) i - 


V Z 2 i^ 



(7) 



'in parallel relay-eavesdropper channel if all subchannels are degraded, the entire relay-eavesdropper channel may not 
necessarily be degraded. 
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where i is the time index, {Zu i}, {Zii} and {Z 2 i,j} are noise processes, independent and 
identically distributed (i.i.d) with the components being zero mean Gaussian random variables 
with variances af t , of and o\\ respectively, for / = 1, . . . ,L. We assume that the source and 
relay know the noise variances present at the receivers. For the subchannel /, Xuj and X 2 i^ are 
inputs from the source and relay nodes respectively. The parameter pu indicates the ratio of the 
R-D link signal-to-noise (SNR) to the S-D link SNR and p 2 i indicates the ratio of the R-E link 
SNR to the S-E link SNR for subchannel / respectively. The source and relay input sequences 
are subject to separate power constraints P 1 and P 2 , i.e., 

- L n 

-££*«,]< a, (8) 

i=i i=i 

L n 

-££E[X 2 2 M ] <P 2 . (9) 

1=1 8=1 



B. Lower Bound on the Perfect Secrecy Rate 

For the parallel Gaussian relay-eavesdropper channel ©, we apply Theorem [2] to obtain a 
lower bound on the perfect secrecy rateE 

Corollary 1: For the parallel Gaussian relay-eavesdropper channel ©, a lower bound on the 
perfect secrecy rate is given by 




Proof: The achievability follows by applying Theorem [2] with the choice Ui := constant, 

Vu := X u , V 2l := X 2U X u := X u + \[^X 2U a, := 1 - a h X xl ~ JV(0, a,P u ) independent 



2 The results established for the DM case can be readily extended to memoryless channels with discrete time and continuous 
alphabets using standard techniques [191 Chapter 7]. 



TO APPEAR IN IEEE TRANSACTIONS ON INFORMATION FORENSICS AND SECURITY 



10 



of X 2 i ~ Af(0,P 2 i), where a t e [0, 1] for / e A; and X u ~ Af(0,P u ) independent of X 2l ~ 
•A/"(0, i^z) f° r I ^ -4 C - Straightforward algebra which is omitted for brevity gives (flOl) . □ 

The parameters Pi; and P 2 z indicate the source and relay power allocated for transmission over the 
l-th subchannel. In (flOl) . after some straightforward algebra, the contribution to the equivocation 
of information sent through NF (set A c in Theorem EJ can be condensed by observing that 
we only need to consider mm{J2 ieA c I(X 2 u Y 2l ), J2 leA c I(X 2l ; Y t )} = J2ieA- K X 2h Y 2i ), to get 
higher secrecy rate. A simplified expression for R l ° w is given by 



R 1 ° w = max min ■ 

Ef=l fl!<-Pl>Ef=l P2l<P2, 

0<q,<1, for; = 1, ... , |«4| 



E 



Pu + P%P% + 2^aip 2 iPiiP2, 



>E 



J aiPu \ _J Pu + P21P21 + 2y / a l p2lPllP2l 



1 + • 



E 



c (Pu±PuPu\ c I Pu + P21P21 



'21 



• E 



A; + P21P21 



(11) 



In (fTT|) . for each subchannel [,] + appears because achievable secrecy rate is always non-negative. 

Remark 5: The achievable perfect secrecy rate established in Corollary \T\ can be larger than 
the one obtained by coding separately over different parallel subchannels. 

This remark is elucidated by the following example. 

Example: We consider a deterministic parallel relay-eavesdropper channel with two subchan- 
nels, i.e., L := 2, as shown in Fig. [3l For subchannel 1, the link capacities to the relay, legitimate 
receiver and eavesdropper are given by Cyi j0 = 4, C r i,6 = 3 and C e \ = 2 respectively. For 
subchannel 2, the link capacities to the relay, legitimate receiver and eavesdropper are given by 
C r 2,a = 5, C r2 ,b = 7 and C e2 = 3 respectively. For this channel, achievable rate obtained by 
coding across subchannels is given by 



Re — Him < E^(^«,a — C e i) + , E^(^«,& — C e i) + > 



min{4,5}=4. (12) 
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Similarly achievable rate obtained by coding separately over each subchannel is given by 

2 

Re = mm {(Cri,a — C ei ) + , (C r ifi — C ei ) + } 
i=l 

= min{2,l} + min{2,4} = 3 (13) 
which is clearly smaller than (fTU) . This shows the usefulness of coding across subchannels. 



C. Upper Bound on the Perfect Secrecy Rate 

The following theorem provides an upper bound on the secrecy rate for the parallel Gaussian 
relay-eavesdropper channel. 

Theorem 3: For the parallel Gaussian relay-eavesdropper channel ©, an upper bound on the 
secrecy rate is given by 



R7 = 

{K P , 



max V I(X lh X 2l ; Y{) - I(X U , X 2l ; Y 2l ) (14) 
? l etc P[ } l=1> ... )L j=i 

where the maximization is over [Xu,X2i] ~ jV(0, Kp,) with JCp, = < Kp, : Kp,= 



Pu_ </>/ y/Pu^a i < jj < iL / = I L, with the covariance matrices 

i>i VP11P21 P21 \ — ti — 

E[X 1[ltL] X^ L] }, E[X 2 [i,l]XJ [1jL] ] satisfying © and © respectively. 

Proof: The proof follows from the rate-equivocation region established for the DM case 
in Theorem \T\ Taking the first term of minimization in the bound on the equivocation rate in 
Theorem [TJ we get 

L 

R e < max^I^,^;^ I U t ) - I(V lh V 2l ;Y 2l | U t ) (15) 



1=1 



where U x -)> (V lh V 2l ) ->■ (X u , X 2i ) ->• (Y h Y lh Y 2l ), for I = 1, . . . , L. The rest of the proof uses 
elements from related works in and 0. Continuing from ([TBI) , we obtain 



(a) 



1=1 
L 



< J2 HVu, V21; Y) - I(V U , V 2l ; Y 2l ) 
1=1 

L 

< 1(V U , V 2l - Y h Y 2l ) - I(V U , V 2l ; Y 2l) 
1=1 
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L 



(b) 



Y}.I{Xx h X 2l] Y h Y 2l ) - I{X lh X 2l -Y h Y 2l | V lh V 2l )\ 



1=1 



-[I{X lh X 2l ;Y 2l ) - I{X lh X 2l ;Y 2l \ V lh V 2l )} 



L 



= X)[J(Xij, Xu] Y h Y 2l ) - I(X lh X 2l ; Y 2l )] 



1=1 



■[I(X U , X 2l ; Y u Y 2l | V 1U V 2l ) - I(X lh X 2l ; Y 2l \ V lh V 2l )} 



L 



<Y,[I(X u ,X 2l ; Y h Y 2l ) - I(X lh X 2l ; Y 2l )] 



i=i 

L 



= Y,I{XiuX 2l -Y l \Y 2l ) (16) 
i=i 

where (a) follows by noticing that I(Vu, V 2 \\ Y\ \ Ui) — I(Vu,V 2 i;Y 2 i | U{) is maximized by 
Ui := constant and (b) follows from the Markov chain condition (Vu,V 2 i) — > (Xi h X 2i ) — > 
(Y h Y Ul Y 2l ),l = l,...,L. 

We now tighten the upper bound (fl6l) by using an argument previously used in 0, in the 
context of multi-antenna wiretap channel. More specifically, observing that, the original bound 
C[5]) depends on p(y h y 2 i\x u , x 2i ) only through its marginals p{yi\x lh x 2l ) and p{y 2i \x lh x 2i ), the 
upper bound (fl6l) can be further tightened as 

L 

R e < min max Y^I{X lh X 2l -Y( | Y^) (17) 

{p(y'i,y' 2 i\ x n> x 2i)} {p{xu,x2i)} 



1=1 



where the joint conditional p{y[, y 2l \xu, x 2 {) has the same marginals as p{y h y 2 i\xu, x 2 {), i.e., 
p(y\\xu,x 2 i) =p(yi\x lh x 2 i) and p(y' 2l \x lh x 2 i) = p{y 2 i\x lh x 2i ). 

It can be easily shown that the bound in ( TT71) is maximized when the inputs are jointly Gaussian, 
i.e., [X lh X 2l ] ~ AT(0, K P ,), K P/ G K Pl with K Pl = |k P; : K P , = * v^OT] j _ x < 

^)/ < l|, / = 1,...,L with the covariance matrices EfX^i] JT^ x j] and E^p^jX^ x j] 
satisfying fO) and © respectively 0,0. 

Next, using the specified Gaussian inputs, and proceeding as in 0, ll20ll . the evaluation of 
the upper bound (fT7l) minimized over all possible correlations between Y{, Y^, for I — 1, . . . , L 
yields 

L 

R e < max Y.IiXiuXrtYfi-IiXuXrtYv). (18) 
This concludes the proof. □ 
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The computation of the upper bound ([141) is given in Appendix HO 

Remark 6: Viewing our Gaussian model © as a specific MIMO relay-eavesdropper channel 
(i.e., one without interference), one can establish a genie-aided upper bound on the secrecy 
capacity of the model © by using recent results on MEMO wiretap channels H-flU, by upper 
bounding the secrecy rate that can be conveyed by the source and relay to the legitimate receiver 
on the multi-access part of the channel with that of an interference-free MIMO wiretap channel 
with 2L-transmit antenna at the sender, L-receive antenna at the legitimate receiver and L-receive 
antenna at the eavesdropper. However, in contrast to CHJ), the upper bound obtained this way 
does not show any degree of separability. More specifically, using H-flU, one can argue that 
the following is an upper-bound on the secrecy capacity of the model ©, 



for some [X 1[ltL] X 2 [i t L]} ~ A/"(0, K P ), and K P = E[(X 1[ltL] X2[i,L])(X 1[hL] X2[i,L]) T } has diagonal 
entries that satisfies ([8]) and © respectively. 

Because the equivalent MIMO channel is interference-free, the upper bound (fl9l) can be written 
equivalently as 



Now, observe that (1201) does not show any degree of separability as in (fl4)) . basically because 
of the additional conditioning on Kj 4-1 , for I — 1, . . . , L. 

Also, investigating our proof in the Gaussian case, one can see that the RHS of (fT3T) and its 
proof are fundamental. As mentioned in the proof, we could obtain the final form (fT8l) essentially 
because the upper bound (fl"5l) that we established depends on the conditional joint distribution 
p{yi,y2i\xu,x 2 i) only through its marginals. 

Example Application: We consider a parallel relay channel with interference at the eavesdrop- 
per. The received signals at the relay, destination and eavesdropper are given by 



Re < ^(^l[l,L]j^2[l,L]; Y[1,L]) — ^(^l[l,L])^2[l,i];^2[l,i]) 



(19) 



L 



(20) 



i=i 




L 




(21) 



k=l,k^l 



V 

interference 
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This model can represent the equivalent channel model of a MIMO relay-eavesdropper channel 
with the interference at the relay and legitimate receiver avoided through singular-value 
decomposition; as the source can always get some feedback from both the relay and legitimate 
receiver, and the relay from the legitimate receiver, which then transforms the MIMO transmission 
into one on parallel channels among the source, relay and legitimate receiver. The eavesdropper 
however does not feedback information on his channel, and so is subjected to cross-antenna 
interference. Constraining the eavesdropper to treat the cross-antenna interference as independent 
noise, one can obtain an upper bound on the secrecy capacity of the model with constrained 
eavesdropper by direct application of (fi4)) . Straightforward algebra gives 



Then, it is clear that the upper bound (1221) holds also for the model (|2TT) with a non-constrained 
eavesdropper. 

D. Secrecy Capacity in Some Special Cases 

We now study the case in which the S-R links are very noisy, i.e., the relay does not hear the 
source. 

Theorem 4: For the model £7]), if the relay does not hear the source: 
1) An upper bound on the perfect secrecy rate is given by 




-l<1pl<l 

for I = 1, . . . , L 




(22) 




(23) 




2) A lower bound on the perfect secrecy rate is given by 




(24) 
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where the maximization is over {Pu,P2i}, I = l,...,L, such that ^2f =1 Pn < Pi, 
Ef=i^<P 2 and 




Proof: 



Upper Bound. The bound in (12~3l) is established as follows. Our approach borrows elements 
from an upper bounding technique that is used in [18], and can be seen as an extension of it 
to the case of parallel relay-eavesdropper channels. Assume that all links between the relay and 
the destination are noiseless, and the eavesdropper is constrained to treat the relay's signal as 
unknown noise. As mentioned in lfT8l . any upper bound for this model with full relay-destination 
cooperation and constrained eavesdropper also applies to the model of Theorem HI 

Now, for the model with full relay-destination cooperation and constrained eavesdropper, we 
develop an upper bound on the secrecy rate as follows. In this case, the destination can remove 
the effect of the relay transmission (which is independent from the source transmission as the 
relay does not hear the source), and the equivalent channel to the destination can be written as 



The eavesdropper is constrained in the sense that it is restricted not to decode the relay's signals. 
Mathematically, this can be stated as follows. Let Z' 2l be a random variable that has the same 
distribution as X 2 i and Pz< (z) = Px 2l {z), and represents unknown noise at the eavesdropper. 
The channel output at the constrained eavesdropper is given by 



For the constrained eavesdropper the relay's transmission acts as unknown noise, with the worst 
case obtained with Z' 2l being Gaussian, for I = 1, . . . , L. The rest of the proof follows by simply 
observing that the resulting model (with the worst case relay transmission to the eavesdropper 
and full relay-destination cooperation) is, in fact, a parallel Gaussian wiretap channel, the secrecy 
capacity of which is established in 0, i.e., 




(26) 




(27) 



unknown noise 



L 




(28) 



i=i 
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where the maximization is over Xy ~ Af(0,Pu) and X 2l ~ Af(0,P 2 i), I = 1,...,L, with 
£f =1 P H < Pi and £f =1 P 2J <P 2 . 

Finally, straightforward algebra which is omitted for brevity shows that the computation of 
(EHJ gives (f23]). 



Lower Bound. The proof of the lower bound follows by evaluating the equivocation in 
Theorem [2] with a specific choice of the variables. More specifically, evaluating © with the 
choice \ A C \ := L, V u := X 1U V 2l := X 2U with X u ~ Af(Q, P u ) independent of X 2l ~ jV(0, P a ), 
I — 1, . . . , L and such that (1251) is satisfied, we get the rate expression in the RHS of ((24|) . The 
RHS of (l24l then follows by maximization over all {Pu, P21}, I = 1> • • • , L, satisfying (1251) and 
the total power constraints J2t=i Pu — Pi an( ^ J2t=i P21 < Pi- □ 

Remark 7: The upper (1231) and lower ([24]) bounds on the perfect secrecy rate of Theorem 
4 have same expressions but are maximized over different input sets. These bounds coincide 
only when the inputs ({Pi;, P21}) that maximize the upper bound (1231) also satisfy d25l) . For this 
specific case, perfect secrecy is established and is given by 

'PA J Pi 



C7 s = max^C ^ -C - (29) 

^ V *i : / W 2 i + p2iP2iJ 



1=1 

where the maximization is over {Pu, P 2 i}, I = 1, . . . , L, such that Yla=i Pu — P> Yla=i P21 < P2 
and 

IV. Example Application 

In this section we apply the results which we established for the Gaussian memoryless model 
in section III to study a fading relay-eavesdropper channel. 

For a fading relay-eavesdropper channel, the received signals at the relay, legitimate receiver 
and eavesdropper are given by 

= h sr jXij + Zii 
Yi = hsd^Xij + h r djX 2 j + Zi 
Y 2: i = h se: iXi : i + h rei iX 2: i + Z 2j i (31) 
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where i is the time index, h sd)i , h rd)i , h seji , h re>i and h sr>i are the fading gain coefficients associated 
with S-D, R-D, S-E, R-E and S-R links, given by complex Gaussian random variables with zero 
mean and unit variance respectively. The noise processes {^i,?}, {Zi}, {Z 2 ,i\ are zero mean i.i.d 
complex Gaussian random variables with variances erf, a 2 and of respectively. The source and 
relay input sequences are subject to an average power constraint, i.e., ^" =1 E[|| X lti \\ 2 ] < nPi, 
^" =1 E[|| X 2 .i || 2 ] < nP 2 . We define hi := [h sdji h rd)i h se>i h rf , ti h STii ] and assume that perfect 
non-causal channel state information (CSI) is available at all nodes. For a given fading state 
realization hi, the fading relay-eavesdropper channel is a Gaussian relay-eavesdropper channel. 
Therefore, for a given channel state with L fading state realizations, i.e., h = {hi}f =1 , the fading 
relay-eavesdropper channel can be seen as a parallel Gaussian relay-eavesdropper channel with 
L subchannels. The power allocation vectors at the source and relay are denoted by Pi(h) and 
P 2 (h) respectively. The ergodic achievable secrecy rate of the fading relay-eavesdropper channel 
(13TT) . which follows from (fTTI) is given by 



i?i ow = max min <! E 



E[Pi(h)]<Pi, 
E[P 2 (h)_]<P 2 , 
0<a(h)<l 



heA 



2C\ 



\h sd \ 2 Pi(h) + \h rd \ 2 P 2 (h) + 2^/a(h)\h sd \ 2 P 1 (h)\h rd \ 2 P 2 (h) 



2C 



2C 



|/l 8e | 2 fl W + \h re \ 2 P2(h) + 2 y/a (h) \h se (h) \ h re | 2 P 2 (h) 



-heA 



2C 



Ih^P^h) + \h re \ 2 P 2 (h) + 2^a(h)\h se \ 2 P 1 (h)\h re \ 2 P 2 (h) 



min < E 



h£A c 



2C ( \h sd \ 2 P 1 (h) + \h rd \ 2 P 2 (h) \ 2C f |fe se | 2 Pi(fe) + |fe re | 2 P 2 W 



a(h)\h ar \ 2 Pi(h) 



^heA" 



2C \hs d \ 2 Pl(h) \ 2C \h re \ 2 P 2 (h) \ 2C ( \h se \ 2 P 1 {h) + \hre?P2{h) 



(To 



(32) 



The upper bound for the fading relay-eavesdropper channel (|3Tb follows directly from the upper 
bound established for the parallel Gaussian relay-eavesdropper channel (fl4l) . Straightforward 
algebra which is omitted for brevity gives 



= max E h { 2C\ 

K[Pi(h)]<P u I ' 
E[P 2 (fe)]_<P 2 , 
-l<ip(h)<l 



\h sd \ 2 Px(h) + \h rd \ 2 P 2 (h) + 2^(h)^\h sd \ 2 P 1 (h)\h rd \ 2 P 2 (h) 



2( , f |fe se | 2 Pl(/l) + \h re \ 2 P 2 (h) + 2 1 p(h)y / \h se \ 2 P 1 (h)\h re \ 2 P 2 (h) \ 1 
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V. Numerical Results 

In this section we provide numerical examples to illustrate the performance of fading relay- 
eavesdropper channel. We consider a fading relay-eavesdropper channel with L realizations of 
fading state. It is assumed that perfect channel state information is available at all nodes. We can 
consider this channel as a Gaussian relay-eavesdropper channel with L subchannels. Alternatively, 
this model can be seen as an OFDM system with L sub-carriers. We model channel gain between 
node % G {s, r} and j G {r, d, e} as distance dependent Rayleigh fading, that is, h it j = h'^jd^J^ 2 , 
where 7 is the path loss exponent, dij is the distance between the node i and j, and h'^ is 
a complex Gaussian random variable with zero mean and variance one. Each subchannel is 
corrupted by additive white Gaussian noise with zero mean and variance one. Furthermore, 
for each symbol transmission same subchannel is used on S-R and R-D links to make the 
optimization tractable. The objective function for both lower and upper bounds are optimized 
numerically using AMPL with a commercially available solver, for instance SNOPT. 

To illustrate the system performance, we set the source and relay power to 64 Watt each. We 
consider a network geometry in which the source is located at the point (0,0), the relay is located 
at the point (d,0), the destination is located at the point (1,0) and the eavesdropper is located at 
the point (0,1), where d is the distance between the source and the relay. In all numerical results 
we set path loss exponent 71=2 and L := 64. For all numerical examples, secrecy rate is given 
by bits per channel use. For each subchannel the selection of the coding scheme at the relay is 
based on the relative strength of the S-D link w.r.t the S-R link, i.e., we use NF scheme (set 
A c ) when \h s d\ 2 > \h sr \ 2 and DF scheme (set A) when \h s d\ 2 < \h sr \ 2 . Fig. 0] shows the power 
allocation for a fading channel with 64 subchannels where the relay is located at (0.5,0), and 
marker 'x' denotes NF on a particular subchannel while marker '□' denotes DF on a particular 
subchannel. It can be seen from Fig. 0] that, achievable perfect secrecy rate is zero for some 
subchannels. Roughly speaking, this happens when the condition \h r d\ 2 > \h re \ 2 is violated. 

Fig. \5\ compares the average perfect secrecy rate of the lower bound, with optimized power 
allocation and with uniform power allocation, i.e., allocating same power at the source and 
relay for all subchannels in h G A and in h G A c . It can be seen that for separate source and 
relay powers, optimized power allocation scheme outperforms uniform power allocation scheme. 
This fact follows because optimized power allocation scheme maximizes the achievable perfect 
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secrecy rate and hence enhances the system performance. 

Mode selection at the relay by only considering the relative strength of the S-D and the S-R 
link in the lower bound is suboptimal because the achievable secrecy rate (|32|) also depends on 
the gain of other link. We now consider the case in which the relay selects the scheme which 
maximizes the rate for each subchannel. We plot the lower bound with this criteria and compare 
it with the case in which same scheme is used on all subchannels. As a reference we consider 
the case in which there is no relay, i.e., a parallel wiretap channel. Fig. [6] shows the achievable 
average perfect secrecy rate of different schemes. It can be seen that when the relay is close to 
the source, DF scheme on all subchannels gives higher secrecy rate. Similarly when the relay is 
close to the destination, NF scheme on all subchannels offers better rate. The region when the 
relay is between 0.5 < d < 1.2 is of particular interest. In this region the relay selects between 
DF scheme and NF scheme for each subchannel and utilizes the gain from both schemes. It 
is interesting to note that when the relay is close to the destination, use of DF scheme on all 
subchannels does not offer any gain because in this case the relay is unable to decode the source 
codewords and hence the average secrecy rate decreases. The lower bound always perform better 
than the wiretap channel which shows the usefulness of the relay. 

In Fig. |7] we compare the lower bound obtained in Fig. [6l with the upper bound on the secrecy 
capacity for the fading relay-eavesdropper channel. It can be seen that when the relay is close 
to the source, the lower and upper bounds coincide. This is achieved by using DF scheme on 
all subchannels. 

VI. Conclusions 

We studied the problem of secure communication over parallel relay channel. Outer and inner 
bounds on the rate-equivocation are established for the DM case. Developing an outer bound 
on the parallel relay-eavesdropper channel is non-trivial and it does not follow directly from 
the one established in lfT2l . For the Gaussian memoryless case, lower and upper bounds on the 
perfect secrecy rate are established. The computable upper bound for the Gaussian model shows 
some separability over subchannels. In the case in which the relay does not hear the source, 
under some specific conditions the lower and upper bounds coincide and secrecy capacity is 
established. We apply the results established for the Gaussian memoryless model to a more 
practical fading relay-eavesdropper channel. Numerical examples showed that power adjustment 
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among parallel channels results in higher secrecy rate. 

Appendix I 
Proof of Theorem 1 

The proof generalizes the results of Theorem 1 in [fT2ll and uses elements from a similar proof 
in the context of parallel BCC in (3). 
a) We first bound the equivocation rate as follows. 

nR e = H(W | Y% 1>L] ) 

= H(W)-I(W;Y 2 n [lM ) 

= I(W; Yfc L] ) - I(W- Y^ L] ) + H(W | Yfa) 

< I(W-Y { l L] )-I(W-Y^ L] ) + ne n 



i=i 

L n 

E 7 ^ y « I Y i~\ Y [i,i-i]) - I(W; Y w I Y? l[i+1] , Y 2 n [l+hL] ) + ne n 



i=i i=i 

L n 



/] E 7 (^' ^2?[i+l] ' ^2[/+l,i] ) I *7 1 ' ^[l,i-l]) ~~ ^(^2?[i+l]) ^2[«+l,L]> I ^> *7 \ ^["/-l]) 
i=l i=l 



L n 

(6) 



i=l i=l 

+ne n 

L n 



J=l i=l 

— /(F/ , yj"/-!]! ^Mi I ^2™[i+l]5 ^2[i+l,L]) ~~ ^(W> I ^1 > *2?[i+l]' ^2[«+l,L]) + ne " 



L n 

— / Yg | 1 ,^[i ] j[_i],^7[i+l]j^2[i+l,i]) ~ I(W;Y 2 u | YJ* \ yji^-i], ^2?[i+l], ^2[/+l 



LI 



i=l 1=1 

+ne n (34) 

where e n — > as rt — > oo; (a) follows from Fano's inequality; and (b) and (c) follows from 
lemma 7 in 0. 
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We introduce a random variable T uniformly distributed over {1,2, ... ,n} and set, Uu := 

Yr\Y { h^ Y W Y Wy V m := W,Y 2 l +1] ,Y? [l+lM and V 2li := l?" 1 , Y^. We define 
17, = (T, C/ji), Vi/ = (T,V m ),V 2l = (T,V 2li ),X n = X lT ,X 2l = X 2T ,Y X = Y T ,Y U = Y 1T ,Y 2l = 
Y 2T , for I = 1, . . . , L. Note that (U h Vu, V 2h X u , X 2h Y h Y lh Y 2l ) satisfies the following Markov 
chain condition 

Ui ->• (V 1U V 2l ) (X u , X 2l ) (Y h Y lh Y 2l ),for I = 1, . . . , L. 

Thus, we have 

- L n 

^ e — ~ / l / l Y H | Yi 1 , y^l-l], *2™[i+l]' ^2p+l,L]) 



1=1 i=l 

—I(W; Y 2U | Yl 1 , ^{"{-i], ^2?[i+i]> ^2[i+i,L]) + e « 

L n 

— ~ / j / j 1 K VV i X l J r 2Z[i+l]' r 2[/+l,L]) *U I Y l 5 r [l,/-l]' r 2Z[i+l]> r 2[/+l,, 
!=1 i=l 

tyw v i_1 v n v n v n ■ v I v n v n v n \ _i_ ^ 

— J^W, Zj , r [i^-i], X 2 l[i+l], *2[l+l,L]l r 2Zi I , 1 j x Z 2Z [j +1 ], I 2 [l+1,L}S + e " 

(J 1 

n 



TV I(V 1H , V 2H ; Y H | C^k) - I(V lH , V 2H ; Y 2H \ U H ) + e n (35) 
n z — ' ^— ' 

1=1 i=l 
L 

HVii, V 2l ; Yi\Ui)- I(V lh V 2l ; Y 2l | U t ) + e n (36) 



(0 

2=1 

where (d) and (e) follow by using the above definition. 

We can also bound the equivocation rate as follows. We continue from d34l) to get 



^ e — n ^ ^ Yh I Y{ 1 , ^"{-i], X2™[i+i]> ^2[2+i,l]) 
2=1 i=i 

— /(W; Y" 2 ,j I ) ^[1,2-1]' *2Z[t+l]> ^2[J+l,i]) + e " 



L n 

-I(W, Yl V , yjy-l], ^[j+l], ^2[i+l,L]; *22i | Yi \ Y[l,2-1]> ^2I[t+l]) ^2[J+l,i]) + e n 
L n 

n ' 

1=1 i=l 

—1{VV, 1; , J: [1,1-1], r 2 l[i+l]j r 2[Z+l,L]i *2H I r « 5 r [1,^-1]' r 2«[j+l]' r 2[i+l,L]J + € n 



~ ^2i[*+l] ' ^2[i+l,L] ) I Yi 1 , ^y-i] , ^2?[i+l] ? ^2p+l, 

«=1 i=l 

ri— 1 -y-n \rn ^rri . y I -yi — 1 yn 

i ) I [l,l-l]i 1 2i[i+l]' - r 2[«+l,L]i J 2ii | *i J 1 [1,1 

- L n 

< — ^ ^2?[i+l] 5 ^2 ? [i+l,L] ) Yin \ Yf 1 , ^[1^-1] , *2Z[i+l] ' ^^[i+l.i] ) 



1 L n 

= TV /(Vik; ii<, Yiu | ^ 2H) - I(Viu, V 2ll ; Y 2U \ U u ) + e n (37) 
n z — ' z — ' 

i=l 8=1 
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(g 

1=1 



1(Vif, Y h Y u | V 2l , Ui) - I(V U , V 2l ; Y 2l | U t ) + e n (38) 



where (/) and (g) follow from the above definition, 
b) We now bound the rate R as follows. 

nR = H(W) 

= I(W;Y { l L] ) + H(W\Y { l L] ) 

< I{W;Y { l L] )+ne n 



ne r 



Y^nw-^riY^) 

1=1 

L n 



1=1 i=l 
L n 



= EE^* i Y r l ' Y [ii-x\)- H ( Y u i WtYr 1 ,^-!]) + 

(■0 



ne r 



1=1 i=l 
L n 

ne r 



1=1 i=l 

U) ^ n 

< EE H <y*) - h v* i w > Y r\ ^-n, Y n^ + ^ 

1=1 i=l 
L n 

— E E ^ X ' ^[l.i-l]' ^2«[»+l]' ^2[Z+l,L]! ^t) + ne n 

(=1 i=l 

L n 

= E E ^ ^) + ne »- (39) 

Hence, we have 

^ L n 

R < - E E 7 (^> ^ *h) + e« 

i=l i=l 

L 

< ^/(^,^;^) + e n (40) 
where (/i) follows from Fano's inequality; (i) and (j) follows from the fact that conditioning 
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reduces entropy. We can also bound the rate R as follows 

nR = H(W) 

= I{W ] Y [ l L] ) + H{W\ Yfa) 

(k) 

< I(W;Y { l L] )+ne n 

L 

= X)W 1 7 1 l i pL,i-i]) + n€ » 

1=1 
L n 



(=1 i=l 
L n 



< E E J ^ ^ ^ i ^ + ne « 

Z=l i=l 
L n 

1=1 i=l 

(0 L n 

< E ^ i ^ ^w-u) - H{ y^ Y » i ^ ^ Y iii-^ y ^+n> + ^ 
i=i i=i 

L n 

- E E ^2J[i+l]' ^2[J+l,i]! ^1£*' I X > ^[M-l]) + n6r 



(=1 i=l 
L n 



E r ^ *« i + ne - (41) 



1=1 i=l 

Hence, we have 

L n 



R ^ - E E i{ y^ Y ^ Y « i + e « 

i=i i=i 

L 

< J2HVu;Y u Y ll \V 2l )+e n (42) 
z=i 

where (k) follows from Fano's inequality; and (/) follows from the fact that conditioning reduces 
the entropy. 

Therefore an outer bound on the achievable rate equivocation region is given by the following 
set: 

[J R e ) that satisfy (EH), (EH), ©J (43) 
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where the union is over all probability distributions p(u [1>L] , v 1[1>L] , f 2 [i,£], x 2[i,l], V[i,l], Vi[i,l], 

U2\i,l})- Finally we note that the terms in (I36l).(l38l).(l40l). and (l42l) depend on the 

probability distribution p(u [1>L] , v 1[1>L] , v 2[1>L] , x 1[1>L] , x 2 [i,l], y[i, L ], 2/ 2 [i,£]) only through 

p(ui,Vii,v 2 i,xii,X2i,yi,yii,y2i)- Hence, there is no loss of optimality to consider only those 
distributions that have the form 

L 

\p{u h v lh v 2 i)p{x lh x 2 i | ui,v lh V2i)p(yi,yii,y2i \ x lh x 2 i). (44) 
1=1 

This completes the proof of Theorem 1. 

Appendix II 

We compute the upper bound on secrecy rate for the parallel Gaussian relay-eavesdropper 
channel as follows. 

L 

max V I(X lh X 2l ; Y t ) - I(X lh X 2l] Y 2l ) 
{Kp,ec P ,} ;= i 

L 

max YlKYi) - h(Y t | X lh X 2l ) - h{Y 2l ) + h(Y 2l \ X u , X 2l )\ 

{Kp ; 6JCp ( } ;= i,... ;i ^ 
L 

max V[/i(y,)-M^i)-M^a) + M^a)]. (45) 
{K P; etp ( }i=i,...,L 

The first term in (|45l ) is computed as follows. 

h(Y t ) = h(X u + yfpTiX 2l + Z,) 

= i log(27re)(P li + Pll P 2l + 2il> iy /p ll P ll P 2l + of). (46) 
Similarly the second, third and fourth term in (|45l) are computed as follows. 

h(Z l ) = Uo g 2TTe(a?) (47) 



(48) 



fr(Y a ) = i log(2 7 re)(P li + p 2i P 2/ + 2^p 2l P n P 2l + 4) 

/i(Z 2{ ) = ilog27re(4). (49) 
Using <|46])-(|49]) in <@3]) gives 

DUP V- 1 , A , p " + Pii^; + 2^i^TWi\ 1 (, Pu+ P21P21 + 2^VP2iPnP 2i \ 
i? e p = max > - log H 5 — - log lH 5 



j:r=iP2i<P2, 
-i<i>i<i 

foil = 1, . . . ,L 



(50) 
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Fig. 2. The parallel Gaussian relay-eavesdropper channel. 
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Fig. 3. An example of a deterministic parallel relay-eavesdropper channel with two subchannels. 
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Fig. 4. Achievable perfect secrecy rate of a parallel relay-eavesdropper channel. 
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Fig. 5. Comparison of achievable perfect secrecy rate of the lower bound with optimized power allocation and with uniform 
power allocation over all subchannels. 
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Fig. 6. Comparison of achievable perfect secrecy rate of some schemes with the lower bound. 
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Fig. 7. Bounds on perfect secrecy rate. 



